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Abstract 

Here we propose a method, based on detrended covariance which we call detrended 
cross-correlation analysis (DXA), to investigate power-law cross-correlations between different 
simultaneously-recorded time series in the presence of non-stationarity. We illustrate the method 
by selected examples from physics, physiology, and finance. 
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There are a number of situations where different signals exhibit cross-correlation. In 
seismology, the degree of cross-correlation among noise signals taken at different antennas 
of detector arrays is an alert signalling earthquakes and volcanic eruptions In finance, 
risk is estimated on the basis of cross-correlation matrices for different assets and investment 
portfolio |2j. In nanodevices for quantum information processing, electronic entanglement 
motivates the computation of current noise cross-correlations, to see whether the sign of this 
signal would be reversed compared to the standard devices |3j. 

Consider two time series {y^} and {y'A-, where i = 1,2, ..,N. Each series can be rep- 
resented as a random walk of k steps, and we can define Rk = y\ + yi + ... + yk and 
R' k = y[ + y' 2 + ... + y' k , where k < N. Series {y^} has a mean // = y~i = YaLiVi 
and a variance a 2 = (y$ — /i) 2 , while series {?/•} has a mean \J = y\ and a variance a' 2 = 
(y'i ~ t 1 ') 2 - We assume that the auto-correlation functions A(n) = (y^ — /x) (yk+ n — f i )/ (j2 
and A'(ri) = (y' k — ^'){y'k+ n ~ t 1 ')/ 17 ' 2 scale as power laws A(n) ~ n~ 7 and A'(n) ~ n -7 ', 
with < 7,7' < 1. We further assume that for the cross-correlation function X(n) = 
(yk — fi)(y'k+ n ~ t L ')/( cro ' f ) between the time series {yi} and {y^} 

X(n)~n" 7x , (1) 

with < 7x < 1. However, this definition assumes stationarity of both time series, and one 
can question its applicability to real-world data typically characterized by a high degree of 
non- stationarity. 

Currently there is no method to quantify the cross-correlations exponent 7 X between 

two correlated time series in the presence of non-stationarity [J]. Here we propose such a 

method, and we illustrate the method by selected examples from physics, physiology and 

finance. To this end, we calculate the covariance: 

n-l 

{R n -TQ{R' n -Rl) =nX{0) + 2Y,[nX{k)-kX{k)\, (2) 

fe=i 

where X(0) = (yk — n)(y'k — fi')/cro J . The sums of Eq. (J2]) can be approximated by inte- 
grals: EkZlX(k) ss Efc=i^" 7x ~ S?dx ar 7x oc n 1 -^ and Y. n k Z\k X(k) S3 ELi ^ 1_7x ~ 
Ji dx x 1_7x oc n 2 ~ 7x . Asymptotically Eq. (jSJ) scales as 

(Rn-TQ^-RZ) ~ n 2A , (3) 

where the scaling exponents A and 7 X - - respectively related to the covariance and the 
cross-correlation function — are not independent, since A = 1 — 0.57 x js|. For {y^} = {y^}, 
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the covariance of Eqs. (2)- (3) becomes the variance that for n » 1 scales as n 2H , so A = H, 
where H is the Hurst exponent. 

In order to quantify long-range cross-correlations when non-stationarities are present, 
we propose a modification of the above covariance analysis which we call detrended cross- 
correlation analysis (DXA). We consider two long-range cross-correlated time series {t/i} and 
{y'j} of equal length N, compute two integrated signals Rk = YJi=i V% an d R'k = Yh=i Vii where 
k = 1, .., N. We divide the entire time series into N — n overlapping boxes, each containing 
n + 1 values. For both time series, in each box that starts at i and ends at % + n, we define 
the "local trend" to be the ordinate of a linear least-squares fit. We define the "detrended 
walk" as the difference between the original walk and the local trend. Next we calculate the 
covariance of the residuals in each box f^xAi 71 ^) = V( n — -0 S/SK-Rfc — Rk,i){R'k ~ R'ki)- 
Finally, we calculate the detrended covariance by summing over all overlapping N — n boxes 
of size n, 

N-n 

F 2 DXA (n)= £ f DX A(n,i). (4) 

i=l 

When only one random walk is analyzed (Rk = R' k ), the detrended covariance F^ XA (n) 
reduces to the detrended variance Fl) FA (n) used in the DFA method [6(. 

In order to test the utility of the proposed DXA method, power-law auto-correlated 
time series yi and y\ are generated by using a stationary linear "ARFIMA" process 
Vi = YljLi a j(p)Vi-j + Vi where < p < 0.5 is a free parameter, aj(p) are weights defined 
by aj(p) = T(j — p)/[r(— p)T(l r(j) denotes the Gamma function, and rji denotes an 

independent and identically distributed (i.i.d) Gaussian variable. The parameter p is related 
to the Hurst exponent, H = 0.5 + p [7\. We generate two time series: {yi} with p = 0.1 
and {y'j} with d = 0.4. Since both {y{\ and {y^} are generated with the same error term 
rji, X(n) ^ |S|. In Fig. 1(a) we show that each time series exhibits the power- law auto- 
correlations expected for ARFIMA, and that the root mean square (rms) of the detrended 
covariance vs. n also follows approximately a power law, consistent with the fact that 
and {y'j} are power-law cross-correlated. For different pairs of power-law auto-correlated 
time series {yi} and {y^}, characterized by Hurst exponents H and H', we find the time 
series are also power-law cross-correlated, where the exponent A of Eq. (3) is approximately 
equal to the average of the Hurst exponents: A ~ (H + H')/2. 

The power-law cross-correlations between {y{\ and {?/•} may exist only if A(n) ~ n~ 7 for 
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both processes. We generate two time series by using two un-coupled ARFIMA processes: 
{yi} with p = 0.1 and {?/•} with p' = 0.4. In Fig. 1(b) we find that, even though both {yi} 
and {y!j} are power-law auto-correlated, the detrended covariance vs. n of Eq. (4) fluctuates 
around zero which indicates that no power-law cross-correlations are present. The same 
result we show for uncoupled ARFIMA processes {y^ and {?/•} defined by p = 0.2 and 
p' = 0.3, respectively. Generally, if the detrended covariance vs. n oscillates around zero, 
there are no power- law cross-correlations with an unique exponent, but either no cross- 
correlations or only short-range cross-correlations exist between {y{\ and {y[}- 

To further exemplify the potential utility of the DXA method for analyzing real-world 
data, we study two time series, both of which can be considered as two outputs of a complex 



system: the air humidity and the air temperature [10] . We analyze absolute values of the 
successive differences of air humidity (denoted by {|?/i|}) and air temperature (denoted by 
{bill) [ see Fig- 2(a)]. Fig. 2(b) shows that each of two time series {\yi\} and {l^'l} exhibits 
power-law autocorrelations with similar scaling exponents. Fig. 2(b) also shows that cross- 
correlations between {\yi\} and {\y'i\} exist and can be fit a power law n x with exponent 
A = 0.75, practically equal to the exponent calculated for the temperature differences. We 
also analyze time series {yi} and {y^}, and we find that the DFA and DXA analyses exhibit 
the correlated behavior, where F^ XA (n) is negative for every n 11 1. 

As a second example of real-world data, we analyze the Sleep Heart Health Study (SHHS) 
database which is designed to clarify the relationship between sleep disordered breathing and 
cardiovascular disease fl^ . For a single patient 12, 13]. we analyze correlations behavior 
of five variables: two electroencephalography (EEG) variables, where EEG is the neuro- 
physiological measurement of the electrical activity of the brain recorded by electrodes com- 
monly placed on the scalp; heart rate (HR) describing the frequency of the cardiac cycle, 
derived from the electrocardiogram (ECG) which records the electrical activity of the heart 
over time; and for both left and right eye the electrooculograms, obtained by measuring the 
resting potential of the retina. 

We analyze the time series of two EEG variables simultaneously recorded every second, 
and find that each of them is short-range auto-correlated. The absolute values of two EEG 
variables we show in Fig. 3(a). The DFA curves in Fig. 3(b) show that each time series of the 
magnitudes exhibits power-law auto-correlated behavior, indicating that a large increment 
is more likely to be followed by a large increment. Fig. 3(b) also shows that, besides auto- 
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correlations, the time series of magnitudes exhibit power-law cross-correlations indicating 
that a large increment in one variable is more likely to be followed by large increment in 
the other variable. We also find that power-law magnitude cross-correlations exist between 
the two time series of magnitudes of two EOG variables. We also find non-vanishing cross- 
correlations between ECG time series and C3/A2 time series. In cross-sectional studies 
where many different physiological time series are recorded, an analysis based on the DXA 
method should add diagnostic power to existing clinical methods employed to discriminate 
healthy from pathological behavior. 

As a third example we analyze the daily closing values of the Dow Jones and the Nasdaq 
financial indices together with their corresponding trading volumes (the number of shares 
traded each day). For both price and trading volumes, we analyze the time series of absolute 
values of the differences of logarithms for successive days. In Figs. 4(a) and 4(c) we show 
their integrated signals I(n) = Z)" =1 (|j/j| — \Vi\)- In Figs. 4(b) and 4(d) we see that each of 
the four time series is power-law auto- correlated, and we also see both for absolute values 
of price changes ( "volatility" ) and trading volume that the time series for Nasdaq and Dow 
Jones indices are power- law cross-correlated. 

We obtain similar results by analyzing cross-correlations between Microsoft and IBM 
stock prices and trading volumes, consistent with the interesting possibility that the above 
results may hold not only for indices, but also for individual companies. This result is 
especially interesting during volatile periods. Long-range cross-correlations between two 
stocks imply that each stock separately has long memory of its own previous values, and 
also long memory of previous values of the other stock. 

Note that it is always possible that cross-correlations between two time series are only 
apparent and exist only due to the presence of long-range auto-correlations in separate time 
series. To test that cross-correlations in magnitudes between two different time series are 
genuine, we generate an artificial time series {xi} that is strongly auto-correlated in the 
magnitudes We find that there are no cross-correlations with an unique power-law 

exponent between an d any of the empirical time series analyzed in the paper. 
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FIG. 1: (a) Rms of detrended variance Fp,FA(n) and detrended covariance, Fr,xA{n), where n is 
a scale. For two time series generated by two ARFIMA processes: {yi} with p = 0.1 and {y[} 
with pf = 0.4 we show the DFA curves Fjj fa (n) for both {yi} and {y'j}, which can be fitted by 
power laws Fdfa ~ n H . Cross-correlations are generated since we choose the error term to be 
equal for both time series: rji = rj^ where rji corresponds to {yi} and rj'i corresponds to {y^}- When 
cross-correlations are present, the same weights are responsible for power-law cross-correlations 
between {y^ and {y^}. For n » 1 we find F^xA(n) ~ n x [see Eq. (4)], where A = 0.73. 
This example illustrates the relation: A ~ (H + H')/2. If we choose the error terms 77- = — rji, 
then Fp XA (n) becomes negative for every n. For that case the cross-correlation function X(ri) 
becomes also negative, (b) Detrended covariance F^ XA (n) of Eq. (4). We generate two pairs of 
two ARFIMA processes, where for each pair the time series are power-law auto-correlated, but 
not cross-correlated, since each ARFIMA is generated by its own error term. The fluctuations, 
both positive and negative, indicate that two time series are not power-law cross-correlated with 
an unique exponent, but either short-range cross-correlated or not at all cross-correlated. 
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FIG. 2: Power-law auto-correlations and cross-correlations in successive differences of air humidity 
{yi} and air temperature {y'i}, recorded each 10 minutes, (a) Time series of their absolute values, 
{\yi\} and {ly^l}. We find that both time series shows sudden bursts of large changes, (b) The rms 
of detrended variance FnFA(n) together with detrended covariance F^xAin). We find that DFA 
curves of {\yi\} and {|?4|} and DXA curve are very similar, and can be approximated with power 
laws Fdfa(ji) ~ n H with scaling exponents H = 0.72 and H' = 0.76, and Fdxa{u) ~ n x with 
A = 0.75. 
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FIG. 3: Power-law auto-correlations and cross-correlations between two different EEG time series: 
CA/Al {{yi}) and C3/A2) ({y'j}) with 32,000 data points, recorded every second, (a) Time series of 
their absolute values, and (b) Rms of detrended variance -Fdfa(^) together with de- 

trended covariance -Fdxa("-)- We find that -Fdfa("-) of {\yi\} and {|^|} and DXA curve, -Fdxa("-) 
are very similar, and can be approximated with power laws n a (n x ) with scaling exponents a = 0.83 
and a = 0.84, A = 0.84, respectively. 
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FIG. 4: Long-range auto-correlations and cross-correlations in absolute values of price changes and 
trading volume for both Dow Jones Index and Nasdaq Index, recorded daily, in the period from 
July 1993 to November 2003. (a) Integrated profiles I(n) = X)£=i(l?/i| — \Ui\) of the time series of 
absolute values {\yi\} and {|y-|} of logarithmic changes in price for Dow Jones and Nasdaq, and 
(c) the corresponding absolute values {\zi\} and {|^|} for trading volume changes, (b) For price, 
rms of the detrended variance Fr,FA(n) curves for each {|^|} and {|z-|}, and also the rms of the 
detrended covariance, Fr>xA{ri). (d) Three curves represents the same measures but for "volume" 
(absolute of trading volume changes) . For all DFA and DXA curves we find both power-law auto- 
correlations and power-law cross-correlations. Power-law cross-correlations between Nasdaq and 
Dow Jones indices imply that current price changes of Nasdaq depend on its previous changes but 
also on previous price changes of Dow Jones. For trading volumes we also analyze time series {zi} 
and {z'j}, and we find that the DFA and DXA analyses show the anti-correlated behavior with 
DFA exponents a = 0.07, a = 0.11, and DXA exponent A = 0.04. 
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